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We report on the electronic ground state of a layered perovskite vanadium oxide Sr 2 V 04 studied 
by the combined use of synchrotron radiation x-ray diffraction (SR-XRD) and muon spin rota¬ 
tion/relaxation (/iSR) techniques, where /xSR measurements were extended down to 30 mK. We 
found an intermediate orthorhombic phase between Tc 2 ~ 130 K and Td ~ 100 K, whereas a 
tetragonal phase appears for T > Tc 2 and T < Td. The absence of long-range magnetic order was 
confirmed by fiSK at the reentrant tetragonal phase below Td, where the relative enhancement in 
the c-axis length versus that of the a-axis length was observed. However, no clear indication of the 
lowering of the tetragonal lattice symmetry with super lattice modulation, which is expected in the 
orbital order state with superstructure of dyz and dzx orbitals, was observed by SR-XRD below Td. 
Instead, it was inferred from /xSR that a magnetic state developed below Tco ~ 10 K, which was 
characterized by the highly inhomogeneous and fluctuating local magnetic fields down to 30 mK. 

We argue that the anomalous magnetic ground state below Tco originates from the coexistence of 
ferromagnetic and antiferromagnetic correlations. 

PACS numbers: 75.25.-j, 75.47.Lx, 76.75.-hi 

I. INTRODUCTION 

The strong spin-orbit interaction is one of the current 
topics in condensed matter physics, as it accompanies 
fascinating phenomena such as the multipole order in /- 
electron systems^ and the Jeff = 1/2 Mott transition in 
5J electron systems^ In this context, the layered per¬ 
ovskite vanadium oxide Sr 2 V 04 is attracting renewed in¬ 
terest as a candidate compound for a Mott insulator 
showing novel spin-orbital magnetism. It has a tetrago¬ 
nal K 2 NiF 4 structure with the space group /4/mmm at 
room temperatureThe magnetic vanadium ion (V^+, 

3d^) is octahedrally coordinated with six ligand oxygen 
ions, where the octahedron is elongated along the c axis. 

The anisotropy of the octahedron leads to a tetragonal 
crystal field that partially lifts the three-fold degeneracy 
of the t 2 g orbitals by splitting them into dxy and dyz/dzx 
orbitals with d^y having a higher energy. Thus, the single 
d electron would occupy the unquenched dyz/dzx orbitals 
in case of a negligible spin-orbit interaction!^ 

However, it has been pointed out that the spin-orbit 
coupling may be large enough to split the dyz/dzx state 
into two Kramers doublets!^ The Kramers doublets are 
mapped to isospin up and down states, which consist of 
the electron spin and orbital angular momentum. It was 
also suggested that a magnetic octupole order (an anti- 
ferroic order of the isospins) would be the ground state of 
Sr 2 V 04 . Furthermore, an alternating scenario of a spin- 


orbital order with almost muted magnetic moments has 
been proposed for the ground state of this compound!^ 

Despite the efforts of several groups to synthesize high- 
quality powdered samples in search for signatures of novel 
spin-orbital magnetism, the varying electronic properties 
reported in earlier studies seem to call for precautions 
in assessing the results in relation to the sample qual¬ 
ity. The magnetic susceptibility (y) and heat capacity 
measurements provided a common ground that the com¬ 
pound successively undergoes three phase transitions at 
Tco - 10 K, Tci - 100 K, and Tc 2 - 130 A re¬ 

port on an in-house x-ray diffraction (XRD) experiment 
showed the splittings of several reflection peaks between 
Tci and Tc 2 ^ It was also reported that a rapid enhance¬ 
ment in the c/a ratio (where c and a are the lengths of the 
unit cell along the respective crystalline axes) occurred 
below Tcii^i^ 

Regarding the magnetic properties, a sudden decrease 
in X was observed at Td, although the magnitude of 
the reduction varied from sample to sample Be¬ 
low Tco, a clear hysteresis was observed in magnetization 
measurements, where the magnetism was confirmed to 
be a bulk property!^ d^d^ The magnetic hysteresis im¬ 
plies the appearance of weak ferromagnetism below Tco- 
Although no experimental evidence for the magnetic oc¬ 
tupole order has been reported so far, recent inelastic 
neutron scattering experiments suggest that the spin- 
orbit coupling is strong enough to induce excitation be- 
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tween the crystal field levelsFurther experimental in¬ 
vestigation is required to clarify the microscopic details 
of magnetism in Sr 2 V 04 . 

A previous study on muon spin rotation/relaxation 
(/iSR) reported a spontaneous internal field at the muon 
site with a broad distribution and residual magnetic fluc¬ 
tuation in the temperature (T) range of 1.8 K < T < 
These results suggest that the spin-orbital order¬ 
ing state, which was proposed by first-principles calcula¬ 
tion studies by Imai et appeared below On 

the basis of the theoretical studies, we can expect a severe 
competition between many types of spin-orbital ordering 
patterns due to the coexistence of ferromagnetic and an¬ 
tiferromagnetic correlations and the associated effect of 
frustration. In such an exotic ground state, it is expected 
that the persistent spin dynamics and short-range corre¬ 
lations are found far below the magnetic ordering tem¬ 
perature. Then, we extended the /iSR experiments down 
to 30 mK. 

In this paper, we report on the electronic properties of 
Sr 2 V 04 inferred from the combined use of synchrotron 
radiation x-ray diffraction (SR-XRD) and /iSR measure¬ 
ments of powdered samples of the same batch. From the 
SR-XRD results, we found an intermediate orthorhombic 
phase in the T range Td <T< Tc 2 , whereas a tetragonal 
phase appears for T > Tc 2 and Td. It is demonstrated by 
/iSR that the anomalies at Td and Tc 2 do not accompany 
long-range magnetic order. Moreover, highly inhomoge¬ 
neous magnetism is observed to develop below Tco- A 
possible model of local held distribution is introduced to 
discuss the observed /iSR lineshape and magnetism be¬ 
hind it over the relevant T region. 


II. EXPERIMENT 

We prepared two batches of polycrystalline samples 
(labeled sample-A and B with a net weight of ^300 
mg and ^30 mg, respectively) by a two step solid- 
state reaction^iS First, the precursor Srio(V 04 ) 6 (OH )2 
was prepared by sintering a stoichiometric mixture of 
Sr(N 03)2 and V 2 O 5 at 550-800 °C for 84 h in air. Then, 
Srio(V 04 ) 6 (OH )2 and SrO were mixed at a molar ra¬ 
tio of 1:2 and sintered at 950 °C for 162 h in a H 2 gas 
flow. /rSR measurements were conducted for sample-A, 
whereas y and SR-XRD measurements were performed 
for both samples, y was measured over a T range be¬ 
tween 2 K and 300 K in the presence of a magnetic held 
of 0.2 T with a superconducting quantum interference 
device (SQUID). 

An SR-XRD experiment was carried out using the 
imaging-plate diffractometer installed on beamline 8 A at 
the Photon Factory at KEK (KEK-PF). The incident 
beam was monochromatized by a Si(III) double crystal. 
X-ray diffraction patterns were measured in the T range 
of 47-305 K by using a He blower with the X-ray wave¬ 
length of 0.077534(3) nm. The lattice constants were ob¬ 
tained by the Le Bail whole-profile fitting to the SR-XRD 


patternsThe Le Bail whole-profile refinements were 
performed using the RIETAN-EP software package^ii 
A conventional /iSR experiment was conducted above 
4.2 K using the DI instrument at J-PARC MUSE and 
down to 30 mK using the DR spectrometer furnished 
with a dilution refrigerator on the MI5 beamline of TRI- 
UME, Canada. The muon spin depolarization function 
G{t) was monitored by measuring the decay-positron 
asymmetry 

A{t) = AoG{t) + Abg, (1) 

where Aq is the initial asymmetry, and Abg is the back¬ 
ground contribution mainly ascribed to the sample holder 
made of silver (in which the muon depolarization is negli¬ 
gible). The /iSR time spectra A{t) was then analyzed by 
curve fitting, where we employed the recursive function 
G{t) derived from the appropriate models for the local 
held distribution (see the next section for more detail). 
We obtained Abg = 0.0221(6) for J-PARC MUSE DI 
and 0.083(3) for TRIUME MI5 DR from the calibration 
measurements. 


III. RESULTS AND DISCUSSION 
A. Magnetic susceptibility 



FIG. I. (Color online) (a) Temperature dependence of the 
magnetic susceptibility for poly crystalline samples of Sr 2 V 04 , 
sample-A (blue) and sample-B (red). The filled and empty 
circles represent data for the zero-field-cooling (ZFC) and 
field-cooling (FC) processes, respectively. The inset shows 
a magnified view near ~I00 K. (b) A magnetization curve at 
2 K for sample-B. 

Figure [TJa) shows the T dependence of y for two poly¬ 
crystalline samples used in the present experiment. We 
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FIG. 2. (Color online) The SR-XRD patterns at 305 K for 
sample-A (blue) and sample-B (red). The green tick marks 
indicate the peak positions of the main phase. The asterisks 
indicate the peaks from the secondary phase. 


observed a sudden decrease in % at Td in both results, 
despite their quantitative differences. We also observed a 
clear difference in the T dependence of magnetization be¬ 
tween the zero-field-cooling (ZFC) and field-cooling (FC) 
processes below Tco- We can see a clear magnetic hys¬ 
teresis in the magnetization curve at 2 K in Fig. Hb), 
suggesting weak ferromagnetism. The saturated ferro¬ 
magnetic moment is estimated to be ^ 3 x IO^^/xb/V. 
The data for sample-A exhibit a small cusp at ^55 K in 
addition to the anomalies at Td and Fco- Although we 
could not identify the origin of the cusp at present, we 
note that the peak temperature is slightly higher than 
the liquid-solid phase transition temperature of oxygen, 
which is usually observed at ^50 K. 

The inset in Fig. [TJa) shows a magnified view near 100 
K, and it clearly shows that sample-B exhibits a larger 
decrease in x at Td than sample-A. It was reported that 
the anomalies at Td and Tco were gradually suppressed 
as the oxygen deficiency S increased for the formula unit 
(Sr 2 V 04 -( 5 )^^ This suggests that sample-A is slightly oxy¬ 
gen off-stoichiometric, although both batches of samples 
were prepared by the same synthetic method. The oxy¬ 
gen off-stoichiometry may be introduced by the prepa¬ 
ration with complex precursor materials and a hydrogen 
treatment. These facts indicate that we need greater con¬ 
trol of the synthesis conditions for more precise control 
of the oxygen content. 
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FIG. 3. (Golor online) Typical examples of synchrotron x- 
ray diffraction data measured at (a) 305 K, (b) 114 K, and 
(c) 66 K. The red lines are the fitted results of the Le Bail 
whole-profile refinement. The green tick marks represent the 
positions of the calculated reflections. The blue lines denote 
the difference between the observed and calculated profiles. 



FIG. 4. (Golor online) Temperature variation of the selected 
peaks in the synchrotron radiation x-ray diffraction patterns. 


B. Synchrotron X-ray diffraction 

The structural properties of the present samples were 
investigated using the SR-XRD technique. Figure [2] 
shows the SR-XRD patterns at 305 K for sample-A and 
B. It has been confirmed that the samples were mostly 
single-phase Sr 2 V 04 with small amount of a secondary 
phase at 20 ^ 15°. The secondary phase was identified 
as Srio(V 04 ) 6 (OH )2 used as the precursor during syn¬ 
thesis. Hereafter, we only show the SR-XRD results for 


sample-B because there are very few differences between 
the two samples. 

Figure [3] shows the diffraction patterns measured at 
(a) 305 K, (b) 114 K, and (c) 66 K with the results of 
the Le Bail whole-profile fitting method. The SR-XRD 
pattern at 305 K was indexed by a tetragonal structure 
with the /4/mmm space group. The lattice parameters 
obtained were a = 0.38402(4) nm and c = 1.25581(8) nm 
for sample-A and a = 0.38393(6) nm and c = 1.25482(8) 
nm for sample-B at 305 K. 

It has been confirmed by the SR-XRD experiment of 
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FIG. 5. (Color online) (a) Temperature dependence of the 
half width at half maximum (HWHM) for selected peaks, (b) 
Temperature dependence of lattice constants. Solid lines are 
guide for the eyes. 


the present sample that a structural transition occurs at 
Tc 2 and Tci, as shown in previous studiesj^i^i^ FigurelHa) 
shows the T variation of the half width at half maximum 
(HWHM) for the selected diffraction peaks. Although 
we could not find any additional peaks in the SR-XRD 
patterns, the HWHM of the 200 and 420 peaks increased 
with decreasing temperature for Td < T < Tc 2 - On 
the other hand, the 110, 004, and 002 peaks did not ex¬ 
hibit line broadening within the experimental accuracy. 
These results might be attributed to the existence of an 
orthorhombic phase for Td < T < Tc 2 . The absence of 
broadening in the 110 peak suggests no diagonal distor¬ 
tion along the (110) direction. It also excludes the coexis¬ 
tence of two tetragonal phases, which had been suggested 
in a previous study4 The enhancement in the HWHM 
for 420(200) peak is explained by the tiny splitting of the 
peak into the 420(200) and 240(020) peaks due to the 
expansion and contraction of the two tetragonal a axes. 

The continuous line broadening and a lambda-like 
anomaly in the T dependence of the specific heat^ 
strongly suggest that the phase transition at Tc 2 is 
second-order. Because no superlattice modulation has 
been discovered below Tc 2 , we can provide the possible 
orthorhombic space group Immm among the subgroups 
of Mimmm. We confirmed a clear difference between 
Mimmm and Immm in the reliability factor (Rwp) of 
the Le Bail analysis for Td < T < Tc 2 , whereas we did 
not observe a large difference between them for T > Tc 2 
and T < Td. For instance, Twp = 9-5% for lAjmmm 
and Twp = 8.6% for Immm at 114 K. The T depen¬ 
dencies of the lattice constants are plotted in Fig. lEb). 
We can see the reduction in the c and a lattice constants 
above Tc 2 , and the a-axis length splits into two below 


Tc 2 . However, the origins of the phase transitions at Tc 2 
are unknown at present. Single-crystal data are highly 
desirable for further understanding. 

Although line broadening was observed for some reflec¬ 
tion peaks for Td < T < Tc 2 as shown above, we could 
not identify clear splitting, broadening, or additional 
peaks in the relevant SR-XRD patterns below Td. The 
HWHMs above Tc 2 and below Td are identical within 
the experimental uncertainty, as shown in Fig. [5](a). Fur¬ 
thermore, we obtained similar values of Twp foi* the Le 
Bail analysis with the M/mmm and Immm structures; 
as shown for T > Tc 2 , Twp = 9.7% for M/mmm and 
Rwp = 9.4% for Immm at 66 K. These results suggest 
that a tetragonal M/mmm phase appears below Td. 

Figures [4](a) and (b) show the 002 and 004 peaks in 
the XRD patterns observed at various values of T. We 
found that the 002 and 004 peaks shifted to lower an¬ 
gles below ^105 K, whereas the 200, 110, and 420 peaks 
moved to higher angles, as shown in Figs. Etc), (d), and 
(f). As shown in Fig. [5] (b), these results correspond 
to a steep enhancement in the c-axis length, whereas 
the a-axis length slightly decreases below Td. This ex¬ 
plains the asymmetric broadening of the 118 and 220 
peaks below Td, as shown in Fig. (U^e) because these two 
peaks move in opposite directions owing to the structural 
change. 

Although there are some anomalies in the structural 
properties, as shown above, we could not find evidence 
of the lowering of the lattice symmetry and superlattice 
modulation below Td, which are expected in the orbital 
order state with the superstructure of the dyz and dzx or¬ 
bitals. This indicates that the degeneracy of the dy^/d^^ 
orbitals is not lifted below Td. These structural prop¬ 
erties below Td are fully consistent with the previous 
results reported by Zhou et alA They suggested that the 
enhancement in the c/a ratio without orthorhombic dis¬ 
tortion might originate from the increased occupancy of 
the dyz/dzx orbitals in place of the dxy orbital. This 
seems reasonable because the crystal field splitting be¬ 
tween the dxy and dyz/dzx orbitals is quite small (^0.08 
eV in the local density approximation calculation) 
Thus, we suggest that the dxy orbital contribution to the 
electronic state is not negligible in the high-T tetrago¬ 
nal phase above Tc 2 owing to the narrow gap between 
the dyz/dzx and dxy states. Then, the electronic state 
changes to the pure dyz/dzx state in the low-T tetrag¬ 
onal phase below Td owing to an enhancement in the 
energy gap, which is caused by the increase in the c/a 
ratio to gain crystal field energy. 


C. Muon spin relaxation 

1. Anomalous paramagnetism (T > Tco) 

In the paramagnetic state of ferromagnetic and anti¬ 
ferromagnetic systems, the muon spin depolarization at 
zero external field has contributions from the static nu- 
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clear dipolar fields and the fluctuating fields from electron 
spins. The nuclear spin contribution is well reproduced 
by the Gaussian Kubo-Toyabe (GKT) function derived 
from the Gaussian distribution of the local field (^loc), 
which is a good approximation for a dense ensemble of 
randomly orientated nuclear magnetic momentsThe 
GKT function is written as 


Ggkt(^) ^ 3 3 


(At)2 


( 2 ) 


where A/y^ = corresponds to the root-mean- 

square for the distribution of the local magnetic field 
(where 7 ^ = 27r x 135.53 MHz/T is the muon gyromag- 
netic ratio). On the other hand, the dynamical electron 
spin contribution is given as exp (—At) with the relax¬ 
ation rate A. Then, the depolarization function G{t) is 
described by the product of the two contributions: 

G{t) = GcKTit) exp(-At). (3) 

Figure [ 6 l^ a) shows the zero-field (ZF)-/iSR spectra at 
300 K and 80 K (i.e., above and below Td), where no 
sign of the emergence of the spontaneous internal mag¬ 
netic field, such as oscillation and/or fast damping, is 
observed below Td. This immediately leads us to the 
conclusion that no long-range magnetic order is present 
below Td, which is in qualitative agreement with the pre¬ 
vious result The temperature dependence of A deduced 
by a curve fitting is plotted in the inset of Fig. [ 6 fa). 

The spectra in Fig. [ 6 l^a) including the one at 18 K 
show that the exponential damping gradually overcomes 
the GKT depolarization below 40 K. Figure [ 6 jb) shows 
/iSR spectra at 18 K measured with ZF and longitudinal 
fields (LFs) of 10 mT and 100 mT. Generally, the expo¬ 
nential damping due to the electron spin fluctuation does 
not depend much on the LF because the fluctuation rate 
of the dynamical field from the ordinary paramagnetic 
moments is much larger than the Larmor frequency of 
the muon spin precession. However, the depolarization 
is quenched by the application of a relatively weak LF, 
which indicates that B\oc is much weaker than the LF 
and that it is quasistatic within the time window of /iSR 
(> 10“^ s). These features allow us to estimate an up¬ 
per limit for the mean hyperfine field at the muon site, 
Tioc — A/T/x ^ 2 mT, and that of the fluctuation rate 

< 0.5 MHz. We also performed curve fitting for the 
LF dependence of the /iSR spectra at 18 K by replac¬ 
ing the exponential term in Eq. (|3]) with the Lorentzian 
Kubo-Toyabe (LKT) function for a finite LF: 

Glkt(^) =1- 

-[H-(-)"]A f (4) 

^0 Jo 

where ujq = is the Larmor frequency of the muon 

spin precession in the presence of an external field Hq, 



FIG. 6. (Color online) (a) Zero-field (ZF)-/iSR time spectra at 
300 K, 80 K, and 18 K. The solid curves are the fitted results 
of Eq. (|3}. The inset shows the temperature dependence of 
the relaxation rate, (b) /xSR spectra at 18 K measured with 
ZF and longitudinal fields (LF = 10 mT and 100 mT) with 
the results of the curve fitting with Eq. (0). 


and jo and ji denote spherical Bessel functions>i^ We 
note that the results of the curve fitting indicate that 
the entire volume of the sample exhibits a Lorentzian- 
like field distribution. Thus, we are led to the conclusion 
that the dilute magnetic moments of d electrons having 
an unusually quasistatic characteristic make a significant 
contribution to the muon depolarization for Tco < T < 
Tci, although their origin is not clear at this stage. 


2. Inhomogeneous magnetism (T < Tco) 

As is observed for the ZE-/iSR spectra in Eig. [7{a), a 
rapid depolarization develops with decreasing tempera¬ 
ture below Tco- We show the ZE- and LE-/iSR spectra 
at 30 mK measured in the presence of various LEs up to 
100 mT in Eig.[3^b). The ZE spectra clearly demonstrate 
the absence of conventional long-range magnetic order at 
the base temperature. The initial fast depolarization ex¬ 
hibits a gradual suppression with increasing LE and is 
nearly quenched at LE = 100 mT. Surprisingly, we also 
observed a finite slope due to the residual magnetic fluc¬ 
tuation in these spectra at a large time, which is almost 
independent of the LE, even at 30 mK. These features 
can be observed from the fast damped 2/3 amplitude 
followed by a slowly relaxing 1/3 component originating 
from the magnetic ordered state with slow persistent spin 
dynamics. 

In a previous /iSR study, Sugiyama et al. analyzed 
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The GbG function is described as one of the general¬ 
ized forms of the GKT function, where A in Eq. m has 
a distribution with the mean value Aq and root-mean- 
square width W. Gchcit) for ZF is written as 

where R is the ratio of the distribution width, and 
R = W/Aq. For = 0, the GbG function is identical 
to the GKT function with a deep minimum after the ini¬ 
tial depolarization. The depth of the minimum decreases 
with increasing R and becomes almost independent of R 
for R > 1. Thus, we imposed the condition 0 < < 1 in 

the following analysis. We note that the GbG function 
does not approach the LKT function, even for i? ^ 1. 

For the curve fitting, we presume that G{t) is primar¬ 
ily described by the product of the GbG function and a 
slow exponential decay due to the persistent spin dynam¬ 
ics in addition to a supplementary term for exponential 
damping: 


FIG. 7. (Color online) (a) ZF-/iSR time spectra in the temper¬ 
ature range between 4.2 K and 10 K. (b) /xSR time spectra 
measured at 30 mK in the presence of various longitudinal 
fields (LFs) up to 100 mT. The inset shows the earlier part of 
the spectrum at ZF. The solid lines are the fitted results by 
using Eq. ®. 


their ZF-/iSR spectra by using a combination of an expo¬ 
nentially relaxing cosine oscillation for the damped 2/3 
component and an exponential function for the 1/3 tail 
signalThey showed that the ratio of the relaxation 
rate of the 1/3 tail component Ataii to the precession 
frequency cjaf is close to unity, Ataii/^AF ^ 1, where 
Ataii/^AF corresponds to the field distribution normalized 
to the mean value of 5ioc- However, their fitting function 
is based on an unique precession signal at a spontaneous 
local field for the conventional magnetic ordered state, 
Ataii/^AF ^ 1 • The observed broad Hioc distribution 
and residual magnetic fluctuation strongly suggest that 
a highly inhomogeneous magnetic state appears in the 
ground state of Sr 2 V 04 . 

The inset of Fig. mb) shows the earlier part of the 
ZF-/iSR spectrum at 30 mK. As mentioned above, the 
GKT and LKT functions are widely attributed to the 
/iSR spectra for random magnetism in dense and diluted 
spin systems. However, the ZF spectrum at 30 mK lacks 
a decrease that usually precedes the asymptotic tail for 
the GKT function, although it exhibits a Gaussian-like 
initial depolarization. We attributed this to a broad dis¬ 
tribution of Hioc, which is not characterized by a sin¬ 
gle root-mean-square width, and adopted the Gaussian- 
broadened Gaussian (GbG) function as a recursive func¬ 
tion for the curve fitting 


G{t) — fGGhG{t)e -h (1 — f)e (6) 

where / is the volume fraction exhibiting inhomogeneous 
magnetism, and A is the relaxation rate reflecting the spin 
dynamics. The second exponential term with the relax¬ 
ation rate Aext is for extrinsic contributions originating 
from the magnetically ordered state region with relatively 
fast persistent dynamics, which makes the above muon 
spin depolarization model unreasonable. The T depen¬ 
dence of / is plotted in the inset of Fig. [HI The inhomoge¬ 
neous magnetism develops over the almost entire volume 
of the sample below 6 K. We obtained R = 0.7(1) be¬ 
low 4.2 K, whereas R is unity at higher temperatures. 
This discrepancy might be caused by the lack of data for 
t ^ 0.1 /is that is masked by the limited time resolu¬ 
tion of the pulsed muon beam (^100 ns). We note that 
the dynamical LKT function and dynamical spin glass 
modelsi^ could not reproduce the spectra over the early 
time range. 

The absence of decrease in the /iSR time spectra 
tells us important aspects of the magnetic ground state. 
Noakes et al reproduced the shallow decrease in the GbG 
function by assuming a local field distribution having ex¬ 
cess low-field sites This implies that the magnetic 
ordered state is accompanied by a large modulation of 
the ordered moment magnitude. They also investigated 
the relationship between R and the correlation length by 
Monte Carlo simulation. On the basis of their results, 
R ^ 0.7 corresponds to a correlation length of ^4 nm. 
Imai et al predicted that the spin and orbital correlation 
functions are reduced as the system size increases 
This result suggests that the correlation length of the 
spin and orbital degrees of freedom are finitely limited. 
Our result strongly supports their predictions. 
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FIG. 8. (Color online) (a) Crystal structure of Sr 2 V 04 with 
the minimum regions of the electrostatic potential. The yel¬ 
low regions define the expected muon site from our Hartree 
potential calculation by VASP code.— This figure is drawn 
by VESTA.— (b) Temperature dependencies of Ao (left axis) 
and the relaxation rate A (right axis). The inset shows the 
temperature dependence of the volume fraction /, showing 
the inhomogeneous magnetism. 


The T dependence of Aq and A are shown in Fig. Mjo). 
Ao behaves like an order parameter below Tco, approach¬ 
ing 10.3(6) at 30 mK. From this value, we can ex¬ 
tract a mean field at the muon site, Aq/t^ = 12.1(7) 
mT. According to our Hartree potential calculation by 
using the Vienna ab-initio simulation package (VASP)^ 
the potential minima for a muon are located around the 
(0.50, 0.24, 0.69) position, as shown in Fig. [8l^a). As¬ 
suming randomly oriented l/ip V d-electron moments, 
we obtained ^52 mT/fiB at the relevant muon 

site. Thus, we roughly estimate the ordered moment at 
the V site to be ^ 0.2/iB- This value is much smaller than 
the dipole moment for the V^+(3d^) ion with quenching 
of the orbital moment, /ieff ^ 1-7/iB for S = 1/2. On 
the other hand, A slightly decreases with decreasing T 
indicating the gradual suppression of the persistent dy¬ 
namics. However, we observed finite residual magnetic 
fluctuation even at 30 mK. Persistent spin dynamics have 
been observed for the magnetic ordered state of several 
geometrically frustrated magnetic systems 

Finally, we briefly discuss the origin of the inhomo¬ 
geneous magnetism. As shown in Fig. [U^a), the mag¬ 
netization data suggest the existence of the oxygen off- 
stoichiometry in sample-A, which is used in the present 
/iSR experiment. However, we could not find any line 
broadenings and additional peaks in the SR-XRD pat¬ 
tern for sample-A, as shown in Fig. [2j This indi¬ 
cates that there are no severe deformation and crystallo¬ 
graphic randomness due to the oxygen off-stoichiometry 
in sample-A. From these results, we can consider that 
the inhomogeneous magnetism might have an intrinsic 
origin, although we would observe a qualitative differ¬ 


ence in the field distribution, correlation length, and re¬ 
laxation rate depending on the magnitude of the oxygen 
off-stoichiometry. Imai et al. also revealed that there 
is severe competition between many spin-orbital order 
patterns due to the coexistence of ferromagnetic and an¬ 
tiferromagnetic interactions and the associated effect of 
frustration in Sr 2 V 04 J^^i^^ In their proposed spin-orbital 
ordering patterns, orthogonal and not strictly orthogo¬ 
nal nearest-neighboring orbital configurations appeared. 
Such orbital configurations induce the coexistence of fer¬ 
romagnetic and antiferromagnetic signs of the exchange 
interactions because of the Goodenough-Kanamori rule. 
We suggest that the observed broad Hioc distribution, 
short-range correlations, and persistent spin dynamics 
are signatures of the severe competition between the pos¬ 
sible spin-orbital ordered patterns in Sr 2 V 04 . 

IV. SUMMARY 

We have conducted SR-XRD and fiSR measurements 
of polycrystalline samples of Sr 2 V 04 . The SR-XRD re¬ 
sults suggest the existence of an intermediate orthorhom¬ 
bic phase between Tc 2 ^ 130 K and Td ^ 100 K, whereas 
a tetragonal phase appears for T > Tc 2 and T < Td. Al¬ 
though a large enhancement in the c/a ratio had been 
confirmed below Td, we could not find clear evidence for 
the superlattice modulation and the lowering of the lat¬ 
tice symmetry from tetragonal /4/mmm, which is easily 
expected for the orbital order state with a superlattice 
structure of dyz and dzx orbitals. These results indicate 
that the degeneracy of the dyz/dzx orbitals is not lifted, 
even below Td. 

The /iSR data revealed that the structural transition 
at Td does not accompany magnetic ordering. On the 
other hand, we extended the /iSR measurement down 
to 30 mK and found a highly inhomogeneous magnetic 
ordered state with slow persistent spin dynamics below 
Tco ^ 10 K. In the magnetic ordered state, the /iSR spec¬ 
tra were well fitted by the GbG function, which is known 
to describe the muon spin depolarization due to the inho¬ 
mogeneous magnetism with short-range correlations and 
a vast moment size distribution. These results suggest 
that the inhomogeneous magnetic state may be a signa¬ 
ture of the severe competition between the many possible 
spin-orbital ordered states predicted by first-principles 
calculation o^^T^^ . 
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